Forces and conductances in a single-molecule bipyridine junction 
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Inspired by recent measurements of forces and conductances of bipyridine nano-j unctions, we have 
performed density functional theory calculations of structure and electron transport in a bipyridine 
molecule attached between gold electrodes for seven different contact geometries. The calculations 
show that both the bonding force and the conductance are sensitive to the surface structure, and 
that both properties are in good agreement with experiment for contact geometries characterized by 
intermediate coordination of the metal atoms corresponding to a stepped surface. The conductance 
is mediated by the lowest unoccupied molecular orbital, which can be illustrated by a quantita- 
tive comparison with a one-level model. Implications for the interpretation of the experimentally 
determined force and conductance distributions are discussed. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 



With the advance of new experimental techniques 
for the fabrication of atomic-scale contacts, it has be- 
come possible in recent years to study electron transport 
through a few or even a single molecule suspended be- 
tween metal electrodes. [J, Q The atomistic details 
of the local contact geometry are crucial in such experi- 
ments, as is evident from the low degree of reproducibil- 
ity of for example the measured I-V curves p| or con- 
ductance traces Since the exact contact geometry is 
difficult to control and characterize, theoretical calcula- 
tions are needed for interpreting and understanding the 
results of the experiments. 

Recently, the conductance of bipyridine molecules in 
a toluene solution trapped between two gold electrodes 
has been measured simultaneously to the force required 
to break the contact, p The statistical data derived from 
~1000 independent measurements show a quantization of 
both the conductance and the rapture force, indicating 
that the contacts consist of single molecules. Such bipyri- 
dine junctions are of interest for two reasons: i) it is one 
of the few systems where fairly reliable measurements 
exist of several properties (forces and conductances) for 
junctions of single aromatic molecules, and ii) pyridine 
contacts represent an alternative to the more commonly 
used thiol anchor groups p. Conductances for a bipyri- 
dine molecule between defect-free Au (111) surfaces have 
recently been calculated by Hou et al. p| for a single 
structure using a cluster model for the description of the 
electrodes. However, in the experiment the surface struc- 
tures are not known, and here we take up the issues of 
stability and the extent to which the force and conduc- 
tance distributions can be related to local structure. 

We have performed plane wave-based density func- 
tional theory (DFT) calculations for the conductance 
and breaking force of a bipyridine-contact. The contact 
between the surface and the molecule is modeled by a va- 
riety of different structures based on the (111) orientation 
of the Au fcc-lattice. By varying the substrate/adsorbent 




FIG. 1: A bipyridine molecule suspended between Au elec- 
trodes. The insets illustrate the variation of the surface con- 
tact geometry, where clockwise from top-right to bottom-left a 
single ad- layer with one to seven Au ad-atoms (white spheres), 
where N c varies from three to nine (see text), on top of a flat 
Au (111) surface are shown. Smaller dark spheres define the 
nitrogen positions. 

distance rigidly we obtain total energy curves, from which 
rapture forces can be derived for the different contact ge- 
ometries. A recently developed scheme for coherent elec- 
tron transport 9] allows us to calculate the conductance 
of the molecule within the same theoretical framework as 
employed for the total energy calculations. Based on our 
data obtained for the different contact geometries we find 
that the force and conductance distributions observed in 
the experiment can only be understood in terms of a high 
structure-selectivity for the contact geometry. 

Fig. n shows an example of the Au/bipyridine/Au 
nano-j unctions and defines the different electrode geome- 
tries investigated in this article. By using our DFT-plane 
wave approach |8| with the PW91 functional for exchange 
and correlation [l0( and an energy cutoff of 340 eV for the 



2 



plane wave expansion, we have established that bipyri- 
dine prefers to bond to the Au (111) surface in an on-top 
position, and therefore all calculations have been per- 
formed for an on-top contact. The supercells consist of 
repeated slabs with seven Au layers in a (111) stacking 
with 3x3 periodicity within the surface plane, where the 
slabs are linked by one bipyridine molecule per unit cell. 
As depicted in Fig. [I] the number of Au atoms in the 
surface plane has been varied from one to seven in such 
a way that the coordination number N c of the Au atom 
directly connected to the molecule varies from three to 
nine. The two structures, where the remaining vacan- 
cies in the surface plane are also filled for iV c =9, are not 
shown in Fig. [3 since our calculations showed that occu- 
pying these atomic positions has no effect on the forces or 
conductances. The geometry of the isolated molecule has 
been relaxed by total energy minimization, but the Au 
atoms have been fixed to their positions in the bulk crys- 
tal structure throughout our study. The torsion angle, 
a, which arises if the two pyridine segments within the 
molecule are rotated relative to each other has been opti- 
mized and found to be around 25° when the molecules are 
adsorbed in the junction. However, since the energy gain 
of this rotation is only around 0.03 eV, we cannot pre- 
dict the effect that a solvent or finite temperatures might 
have on it. We found that the influence of this rotation 
on the calculated forces is negligible. For the conduc- 
tance we have performed calculations for both a=0° and 
25°. An optimization of the Au-N bond length gave val- 
ues between 2.12 and 2.42 A depending on the number 
of atoms in the surface plane but showed no dependency 
on a. 

In the top panel of Fig.[2|the binding energies at the op- 
timal Au-N bond lengths for all the structures are shown. 
The binding energies vary significantly with an almost 
linear dependence on 7V C , and their magnitude indicate 
intermediate-strength bonding stronger than typical van 
der Waals bonding but weaker than covalent bonding be- 
tween open shell systems. Depending on the detailed 
choice of exchange-correlation functional the binding en- 
ergies may shift a few tenths of an eV but the trend with 
metal coordination number is unchanged (see Fig. (2J. 
The same linear dependence of the binding energy on 
N c has been found for the chemisorption of oxygen and 
CO on Au chains and surfaces indicating that the 
nature of the bonding is also covalent in the present case. 
A closer inspection of the surface states shows, that the 
trend in the binding energies correlates with an upward 
shift of the (i-states towards Ep which is expected to en- 
hance the energy gain of a hybridization with the highest 
occupied molecular orbital (HOMO, see inset) located at 
~ -4.6 eV. 13J These findings are in contrast with a recent 
article, in which the adsorption of pyridine on Au (111) 
surfaces was interpreted in terms of Van der Waals inter- 
actions based on more qualitative arguments such as 
the amount of charge buildup in the center of the bonding 



region. 

From the total energy curves obtained by varying the 
Au-N bond lengths, the breaking forces can be calcu- 
lated by taking the derivatives at the inflexion point. 
As can be seen in Fig. El the forces like the energies de- 
pend almost linearly on N c . The rather large variation of 
the calculated breaking forces should be compared with 
the experimental situation where a quite sharp peak at 
a value of 0.8 ± 0.2 nN is found, p This distribution 
is consistent with our results if the experimental situa- 
tion results in the preference for a subset of structures 
in the range N c =A-7. Such a preference could be either 
thermodynamic or kinetic in origin. Energetically bipyri- 
dine clearly favors to bond to metal sites with low metal- 
coordination numbers according to the calculations. On 
the other hand such metal defects may be energetically 
costly to produce and therefore limited in number. Also 
the experimental process in which the nanocontact is cre- 
ated may bring the system out of local equilibrium and 
the resulting structures will then depend on the evolu- 
tion of the process. Finallly it should be noted that the 
calculations do not take potential effects of the solvent 
into account. 

The conductance of the bipyridine junctions is calcu- 
lated using a numerical method for phase-coherent elec- 
tron transport which combines a Green's function for- 
malism with a plane wave-based DFT description of the 
atomic and electronic structure of the entire system. The 
Green's function of the central region, comprising the 
molecule and part of the gold electrodes, is evaluated in 
terms of a basis consisting of maximally localized Wan- 
nier functions which are defined by an appropriate trans- 
formation of the Kohn-Sham eigenstates. Despite their 
optimal localization, the Wannier functions still have a 
fairly long-ranged but small tail, which for the systems 
studied in this article we truncate at a distance of 19.4 A 
from the center of each function. The couplings to the 
semi-infinite gold electrodes are included via self-energy 
matrices which are also represented in the Wannier func- 
tion basis. A detailed discussion of the transport scheme 
can be found in Ref. fej. The central regions used in the 
transport calculation are the same as the supercells ap- 
plied in the total energy calculations. These supercells 
contain 3x3 atoms in the directions perpendicular to 
the transport direction, and we have found that k-point 
sampling in the transverse plane is important in order to 
obtain well converged results for the conductance. Thus 
all calculations presented in the following have been per- 
formed using a 4 x 4 k-point grid in the transverse Bril- 
louin zone. For a systematic study of the role of k-point 
sampling for transport calculations see Ref.GjJ 

In Fig. 03 we have plotted the calculated transmission 
near the Fermi level (which is taken as the zero point) for 
all geometries outlined in Fig. [l] The top panel of Fig. 
shows the conductances for nano-j unctions with N c =3 to 
9 for a=0° and a=25° as black and grey solid lines, re- 
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FIG. 2: Top: Binding energies for a single N-Au bond as 
a function of N c . The inset shows the HOMO. The solid 
(dashed) line marks calculations using the PW91 (RPBE 
pJ| ) XC-functional. The position of the d-band center with 
respect to the Fermi level is shown as a dotted line. Bottom: 
Rapture force as a function of N c . The horizontal dotted line 
indicates the experimental value and the grey shaded area 
marks the region defined by the experimental error bars. 



spectively. The conductance is seen to vary significantly 
with coordination number. However, the conductance 
for a=25° is in good agreement with the experimentally 
observed value of one percent of the conductance quan- 
tum Go = 2e 2 //i.|5| It can be noted that this is the same 
range of N c ~5, where also the forces in Fig. [2| match the 
measurements. But as already discussed in connection 
with the force distributions the actual propensity of the 
different structures in the experiment is not known. 

We now proceed to analyze why the conductance varies 
the way it does with coordination number. The analysis 
is based on the calculations with a = 0°. 

Within our calculational scheme it is possible to define 
molecular orbitals even when the molecule is coupled to 
the metal. The molecular orbitals are constructed from 
Wannier functions located on the molecule |£j . It is there- 
fore also possible to study the transport through a single 
or a subset of the molecular orbitals by simply removing 
the remaining orbitals from the basis set. When we focus 
on the contributions coming from individual molecular 
orbitals, we find that the single peak in the transmission 
function around Ep (bottom panel in Fig. |3J) is exclu- 
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FIG. 3: Top: Conductance in dependence on N c for torsion 
angles a—0° (black line) and 25° (grey line) in units of the 
conductance quantum Go = 2e 2 /h. The dashed line refers to 
a simple model explained in the text. The horizontal dotted 
line shows the experimental value. Bottom: Transmission 
functions for all nano-j unctions (a=0°). The peaks are due 
to the LUMO, shown in the inset, which rises continuously in 
energy as N c increases. 

sively due to the LUMO, whose wavefunction is shown 
in the inset. Since the peak is quite narrow for all geome- 
tries, correlation effects, leading to Coulomb blockade be- 
haviour in other weakly bonded molecular junctions 
could be relevant for its accurate description. However, 
even in the "coherent limit" there is good agreement be- 
tween our transport calculations and the experimental 
results. 

One way of analysing the shape of the T(N C ) curve 
further is by applying a simple one-level model 14] with 
only one molecular state (the LUMO) coupled symmet- 
rically to the leads . Within this model the transmis- 
sion through a level at energy Eq can be expressed as 
T(E) = A 2 /((E - E F ) 2 + A 2 ), where the peak-width is 
denned by A = ttV 2 p° . V is the coupling constant of the 
LUMO to its group-orbital, 0, and p® is the density 
of states of group orbital in the absence of coupling to 
the LUMO evaluated at Eq. By taking the values for 
P®,Eq, and V from the DFT calculation and using them 
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to calculate A, we obtain the dashed line in Fig. 03 which 
thus represents the conductance when the electrons are 
allowed to travel between the leads via the LUMO only. 
The good agreement of this one-level model when com- 
pared with our data further illustrates that the variation 
of the conductance is entirely controlled by the position 
of the transmission peak due to the LUMO and that the 
effect of the interference with other orbitals as well as 
through- vacuum tunneling are negligible. 

But what is governing the LUMO position? We in- 
vestigated the variation of the energetic position of the 
lowest-lying molecular orbital (denoted MOl) with re- 
spect to the metal's Fermi level, in dependence on the 
distance between surface and molecule for all studied 
surface structures. In this analysis it can be seen that 
at the bonding distance MOl and the LUMO depend on 
N c in quantitatively the same way. This pattern cannot 
be explained in terms of vacuum level alignment alone, 
since the variation of the metal's workfunction with 7V C 
is parabolic with a minimum at 7V c =5-6. However, MOl 
is ~ 10 eV below the lowest-lying Au valence states, and 
therefore the exchange or bonding interaction between 
the molecule and the surface do not have a direct influ- 
ence on the energy of this orbital. From that we can con- 
clude that the dependence of MOl as well as the LUMO's 
energy on N c must be exclusively guided by rigid poten- 
tial shifts, where we found that equilibrium charge trans- 
fer dominates at the optimal bonding distance and also 
Coulomb repulsion due to the confinement of electrons 
at the interface plays a role. We will discuss these elec- 
trostatic effects and their structure-dependent influence 
on the aligment of molecular levels in detail in a separate 
publication. 

In summary, we presented a theoretical analysis of the 
forces and conducances of a bipyridine molecular junc- 
tion. Since important parameters of the experiments, 
such as the atomic configuration of the electrode's sur- 
faces are not known, we scanned a wide range of surface 
structures and found that both the forces and conduc- 
tances depend crucially on their geometries, but our cal- 
culations agreed with the experiments in orders of mag- 
nitude for both quantities. The narrow peaks found in 
the statistics of the experiments, however, could not be 
directly reproduced from our calculations, and the statis- 
tics of the measurements site-selectivity of the adsorption 
of the molecules has to play a major role. Our study 
suggests that varying surface structures might be an im- 
portant ingredient for a systematic understanding of the 
bonding and conductance properties of molecular junc- 
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